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We investigate the electronic structure of the transition-metal oxide YVO3 by a hybrid 
first-principles scheme. The density-functional theory with the local-density-approximation 
by using the local muffin-tin orbital basis is applied to derive the whole band structure. The 
electron degrees of freedom far from the Fermi level are eliminated by a downfolding proce- 
dure leaving only the V id t2g Wannicr band as the low-energy degrees of freedom, for which 
a low-energy effective model is constructed. This low-energy effective Hamiltonian is solved 
exactly by the path-integral renormalization group method, ft is shown that the ground 
state has the G-type spin and the C-type orbital ordering in agreement with experimental 
indications. The indirect charge gap is estimated to be around 0.7 eV, which prominently 
improves the previous estimates by other conventional methods. 
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1. Introduction 

Recent intensive studies have revealed that strongly correlated electron systems exhibit 
a variety of distinguished physical properties.^ Basic aspects of the correlation effects have 
been studied mainly by simplified models, such as the Hubbard, and Heisenberg models. For 
example, at least some essence of the Mott transition, which is widely seen in the strongly 
correlated electron systems, is believed to be captured by the single-band Hubbard model. 
However, if we focus on a variety of physical properties in real materials, many challenging 
issues are found beyond the simplified models. Especially in the systems where spin and orbital 
are coupled, rich structure and various phenomena have attracted much interest. When the 
strong Coulomb interaction suppresses charge fluctuations, spin and orbital degrees of freedom 
play a crucial role. In order to understand such realistic and complex systems, computational 
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methods from the first principles offer promising approaches because of their ability to treat 
charge, spin and orbital degrees of freedom on equal footing with full fluctuation effects taken 
into account. 

Computational methods for electronic structure calculations are roughly classified into 
two categories. One is based on the density-functional theory (DFT) supplemented, for ex- 
ample, with the local-density-approximation (LDA).^'^ Since the LDA costs less computation 
time, this scheme has proven its advantage in calculating electronic structure of real ma- 
terials. However, the LDA essentially ignores spatial and dynamical fluctuations which are 
relevant in the strongly correlated systems. In addition, the LDA in general underestimates 
the charge gap. The other approach is based on the wave function scheme. The Hartree-Fock 
(HF) method often gives a good starting point to discuss competition among ordered states, 
while it tends to overestimate the charge gap. The HF approximation itself is a rather crude 
one, where spatial and temporal correlations are ignored. However, one can improve the ac- 
curacy systematically within the wave function scheme by expanding basis set, for instance, 
by introducing configuration interaction or by path-integral renormalization-group (PIRG) 
methods.^ The method along this line tends to need much longer computational time, which 
prevents us from applying it directly to real materials. 

Recently, a hybrid method called DFT-PIRG which combines the DFT with a wave func- 
tion method, namely PIRG has been proposed.^ The DFT is utilized to first obtain the overall 
band structure extending to the energy scale far from the Fermi level. Then the downfold- 
ing procedure eliminates the high-energy degrees of freedom far from the Fermi level by the 
renormalization process. After this downfolding process, an effective model for the low-lying 
electronic structure is obtained, which is solved exactly by a low-energy solver. The PIRG 
method is a promising choice for a low-energy solver because it allows us to treat spatial and 
dynamical fluctuations in a controllable way. It has indeed been applied to Sr2V04^'^ and has 
reproduced the basic experimental properties including the fact that this compound is located 
just on the verge of the Mott transition with a tiny but nonzero gap amplitude (< O.leV). 
In the present work, we study YVO3 as another example of the charge-spin-orbital coupled 
system by the DFT-PIRG method. 

YVO3 belongs to the family of transition-metal oxides with two valence electrons in the 3d 
orbitals {t2g manifold). The lattice structure is an orthorhombically distorted perovskite with 
the space group Pbnm (four vanadium sites in a unit cell) at room temperatures. The GdFeOs- 
type distortion, rotation and tilting of the VOe octahedra are present, where the reduced 
V-O-V angle makes the narrow t2g bands. With lowering the temperature, it undergoes two 
successive phase transitions in both spin and orbital sectors. First, the G-type orbital ordering 
(00) appears at 200K with a structural change to the P2i/a symmetry, where a site with the 
dxy and dyz orbitals occupied and one with the the dxy and d^x are alternatively arranged in 
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three dimensions. The magnetic structure shows the C-type spin ordering (SO) below 116K, 
where spins are ahgned antiferromagnetically in the a-b plane and ferromagnetically along 
the c-axis. With further lowering the temperature, the SO and 00 simultaneously change at 
77K, and the ground-state is the C-type 00 with the G-type SOJ'^ The crystal structure 
recovers the Pbnm symmetry as illustrated in Fig. 1. In the charge sector, YVO3 is a typical 
Mott insulator with a large charge gap (~ leV). This is partly attributed to a large GdFeOs- 
type distortion, which reduces the band width effectively. In addition, coupling to Jahn- Teller 
distortions is important in determining the orbital states. 



Fig. 1. (Color online): Illustration of the Pbnm orthorhombic crystal structure of YVO3. Each V 
atom (large red sphere) is octahedrally surrounded by O atoms (small blue spheres). Other (light 
yellow) spheres represents Y atoms. 

The relation between the SO and the 00 has been studied by the unrestricted HF approx- 
imation.^' ^'^ Although the HF result explains the spin and orbital orders in the ground-state 
consistently with the experiments, the magnitude of the charge gap is largely overestimated (~ 
3.4eV) compared to the experimental value (~ leV). From the first-principles calculations, the 
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local spin density approximation (LSDA) and the generalized gradient approximation (GGA) 
have been applied with the full-potential linearized augmented-plane-wave 

method. ^^'^2 ^j^g 

LSDA results have shown the metallic ground state. On the other hand, the ground state 
is insulating with GGA, but it failed to reproduce the correct SO and 00. In addition, the 
gap amplitude is underestimated (0.009eV). The LDA+C/ method has succeeded in describing 
the correct ground state. ^^'^^ However, aside from the issue of avoiding the double counting 
of the Coulomb interaction generally known in LDA+C/ calculation, this ground state is ob- 
tained by adjusting the parameter U to reproduce the experimental band gap. It is desired 
to estimate the screened Coulomb interaction itself from the first principle. In addition, it is 
also important to go beyond the single-Slater-determinant approximation to examine effects 
of quantum fluctuations accurately, and to estimate spatial and temporal fluctuations from 
the first principles. This is because orbital-spin fluctuations may seriously alter the ground 
state through strong correlations combined with quantum fluctuations. 

The purpose of the present work is to study the ground state of YVO3 from the first- 
principles calculation at a quantitative level. We will discuss how the SO and 00 are repro- 
duced with the DFT-PIRG scheme and estimate the charge gap. The basic procedures are the 
same as the case of Sr2V04. Readers are referred to refs. 5,6 for more detailed description of 
this scheme. Here we note that the target materials are quite different, while the procedures 
are similar. In contrast to Sr2V04, YVO3 is a cubic perovskite with the GdFeOs-type distor- 
tion, where its isotropic three dimensionality or mixing between t2g orbitals may show different 
aspects. The inter-orbital Coulomb interaction and the Hund's rule coupling are considered 
to be more effective for the 3(i^ configuration of YVO3. Furthermore, YVO3 has much larger 
Mott gap (~ leV). To examine the applicability of the DFT-PIRG method, it is desired to 
test whether it works in such completely different systems. 

In Sec. II, we derive the effective low-energy model. In Sec. Ill, the Hartree-Fock solution 
to this model is considered. We discuss results of the PIRG calculation in Sec. IV. Section V 
is devoted to summary. 

2. Constructing Effective Model 

In many cases of the strongly correlated electron systems, the bands close to the Fermi 
level are relatively narrow and well separated from the bands far from the Fermi level. In 
fact, for the transition metal oxides, the Sd bands are mostly responsible for the low-energy 
excitations. Compared to the width of the 2>d bands, other bands are located far from the Fermi 
level. Furthermore, if the crystal-field splitting is strong, it may be sufficient to consider only 
either Cg or t2g bands for the low-energy degrees of freedom. This hierarchy structure in energy 
enables us to treat the system in a hybrid scheme. 

The electron correlation effects and dynamical fluctuations are important near the Fermi 
level, whereas the electronic structure far from the Fermi level is expected to be well repro- 
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duced by LDA. Therefore, in the first step of our hybrid scheme, the whole electronic struc- 
ture is calculated from the density functional approach by employing LDA based on the linear 
muffin-tin orbital (LMTO) method. ^^'^^ The LDA enables us to calculate the whole electronic 
structure within the presently available computer power. Then we eliminate the electronic de- 
grees freedom far from the Fermi level. This is done by the renormalization (downfolding) 
procedure, where the electron correlation effects are properly considered in the calculation of 
the screening and self-energy effects. After the downfolding procedure, the low-energy degrees 
of freedom isolated near the Fermi level are extracted in the effective Hamiltonian, which will 
be solved accurately in the PIRG method. 

Let us explain the whole procedure developed in the literature^' ^ by extending it specifi- 
cally for the present compound YVO3. At first, we compute the LDA band structure from the 
LMTO basis. Since YVO3 has sufficiently large crystal-field splitting (~ O.leV) due to a large 
Jahn- Teller distortion, we focus on three low-energy bands which mainly consist of the 3d t2g 
orbitals. The whole Hilbert space is now spanned by the LDA basis functions for the LMTO 
t2g bands {\d)} and the set of the rest basis {|r)}. Then we eliminate the subspace and 
reduce the Hilbert space to the subspace which only consists of the set The Hamil- 

tonian matrix is now spanned only within this restricted Hilbert subspace, which provides us 
with the low-energy tight-binding effective Hamiltonian defined on the Wannier basis for the 
3d t2g bands. Note that the band structure of the effective tight-binding Hamiltonian defined 
only on the low-energy subspace well reproduces the t2g part of LDA bands computed from 
the whole LMTO basis as is illustrated in Fig. 2. 




r xs YrzuR tz 



Fig. 2. (Color online): Comparison of the band structure of 3d t2g orbitals computed from LMTO 
calculations (solid light blue lines) with the downfolded tight-binding model (dashed-dotted brown 
lines) for YVO3. 

Next, the parameter values for the final effective low-energy Hamiltonian for the 3d t2g 
Wannier bands become renormalized by the interaction between electrons in the subspace 
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{|r)} and those of 3d t2g near the Fermi level in the subspace The renormalization 

appears both in the screening of the interaction parameters and the renormalization of the 
kinetic energy part. These two procedures are briefly outlined in the following paragraphs. 

Here we sketch out the derivation of the effective Coulomb interaction among the 3d t2g 
electrons, which consists of two steps. In the first step, by using the standard constrained-LDA 
(c-LDA) method, we compute the interaction matrix Wri without effects of hybridization be- 
tween 3d and non-3d orbitals. This part takes into account only the screening caused by the 
redistribution of the non-3(i electrons. From this c-LDA procedure, we obtain the screening 
of the Coulomb interaction between two electrons in the t2g Wannier band caused by the 
polarization of electrons at the non-3(i orbitals in the LMTO basis. In the second step, we 
consider the screening of the Coulomb interaction between t2g electrons caused by polariza- 
tion of electrons on the 3d LMTO basis other than the t2g basis (namely caused by the eg 
component and the 3d LMTO atomic orbital component residing in the oxygen 2p Wannier 
band) in the random-phase-approximation (RPA) scheme. By taking Wri obtained from the 
c-LDA scheme as a starting Coulomb interaction, we can write the RPA screening as 

where Pdr denotes polarization of electrons in the 3d LMTO basis but in the non-t2g Wannier 
bands. The total polarization of the 3d bands is given by 

Pd{uj) = Pdr{uj) + Pt2g{u:), (2) 

where Pt2g represents contribution from 3d t2g Wannier bands. Then, we notice that the 
following identity holds, 

^ ^ 1 - Pd{0j)Wrl 1 - Pt2g{0j)Wr{uj) ^ ^ 

It implies that Wr{uj) can be considered to be the effective Coulomb interaction among the 3d 
t2g electrons. Although this partially screened interaction Wr{uj) in general is dynamical, its 
frequency dependence is small within the energy scale of the t2g bandwidth. Thus we can safely 
take the zero-frequency limit Wr{u> = 0) as the effective Coulomb interaction to construct the 
low-energy Hamiltonian, which will be treated exactly by the PIRG method. We note that 
this separability of the screening effects enables the downfolding procedure. The reason why 
we perform two steps in calculating the screening is that the c-LDA is less computer-time 
consuming, whereas it does not take into account the frequency dependence. On the other 
hand, RPA is more time consuming, whereas it is to some extent able to consider dynamical 
fluctuations, which is important near the Fermi level. 

In the downfolding process, the kinetic-energy part is also modified through the self- 
energy S(A;, to). Such a self-energy is mostly momentum independent and its imaginary part is 
negligible at low energy. Thus the self-energy effect mainly contributes to the renormalization 
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factor Z given by, 



Z = [1 - aReS/acu|<^=o] 



-1 



(4) 



which is obtained by the numerical estimates from the w-dependence of the self-energy S(u;). 
This means that the self-energy effect reduces the band width, leading to the kinetic-energy 
part renormalized by Z. The self-energy can be evaluated in the GW approximation, 

and the renormalization factor Z has been calculated for several transition metal oxides. We 
take Z = 0.8, which is a typical value among these materials. For more complete description 
of the downfolding procedure, readers are referred to refs. 5, 6, 17, 18 and 19. 

After eliminating the high-energy part, the effective Hamiltonian is reduced to a three- 
band Hubbard model in three dimensions: 



where TYk and 7i\j denote the kinetic and interaction energy terms, respectively and A is a 
tuning parameter which controls the overall strength of the interaction. The realistic value 
corresponds to A = 1. Due to the lattice distortion, each unit cell has four vanadium sites. 
The definition of the unit cell is shown in Fig. 3. The kinetic part of the Hamiltonian is given 
by a tight-binding model. 



where cli^^{ciima) denotes the creation (annihilation) operator for spin a =], J, at m-th orbital 
of site I in the i-th. unit cell and tiju'mm' is the overlap integral. We note that the orbitals which 
we consider here are "local t2g orbitals" , which is constructed basically from the local crystal- 
field of the VOg octahedral frame. Because of the tilting distortion of the VOg octahedra, 
each orthogonal orbital is a mixture of dxy, dyz, (1^2, dzx and dx2_y2 in the original coordinate. 
However, the t2g components are still dominant due to the large crystal-field splitting. 

We choose a base \l,m) that diagonalizes the on-site matrix elements tmimm'-'^^ The ex- 
pansion coefficients of the unitary transformations over the real harmonics {\xy), \yz), 
\zx), — y^)) in the orthorhombic coordinate frame are given as follows: 



7i = T~L\ + XH\] 



(5) 




(6) 



o- i,j 1,1' 



m,m' 



( V-^y) \ 




(7) 
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Ml 



M2 



/-0.2302 
-0.8913 
\ 0.0605 

/-0.2304 
-0.8912 
\ 0.0603 

/-0.2318 
-0.8909 
\ 0.0601 

/-0.2302 
-0.8913 
y 0.0608 



-0.3619 
0.2435 
0.7922 

0.3617 
-0.2434 
-0.7924 

0.8193 
-0.1044 
0.4719 



-0.3785 
0.0527 
0.2544 

-0.3784 
0.0527 
0.2544 

-0.3784 
0.0533 
0.2544 



0.8195 
-0.1027 
0.4720 

-0.8195 
0.1030 
-0.4718 

-0.3614 
0.2439 
0.7924 



0.0355 \ 
-0.3646 
0.2851 J 

0.0355 \ 
-0.3647 
0.2850 ^ 

-0.0349\ 
0.3648 
-0.2849 J 

-0.0357\ 
0.3646 
-0.2852y 



(8) 



(9) 



(10) 



(11) 

and dyz {dy^ 



-0.8198 -0.3783 0.3614 
Ma = -0.8913 0.1026 0.0527 -0.2437 
-0.4715 0.2546 -0.7924 

where we can see that the orbitals m = 1 and m = 3 mainly correspond to d 
and dzx) orbitals, respectively for / = 1,2(3,4) and the orbitals m = 2 mainly consist of dxy 
The diagonalized matrix has (m, m) elements given by 

/-0.2711 -0.0000 0.0000 \ 
Uiiimm' = 0.0000 -0.2175 0.0000 , (12) 
y 0.0000 -0.0000 -0.1054y 
within each vanadium atom (namely only diagonal elements for i and / are shown here). This 
matrix elements do not depend on i and /. We summarize other off-diagonal matrix elements 
in Appendix . The lattice distortion brings about level splittings among t2g orbitals, which 
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Uii U22 f^33 


U12 Uis U23 


J12 Jl3 J23 


energy (eV) 


3.178 3.235 3.232 


1.961 1.955 1.985 


0.633 0.627 0.633 



Tabic I. The interaction energy obtained by the downfolding procedure. 



is estimated to be 0.17 eV as seen in eq.(12). We take into account up to the third-nearest- 
neighbor transfers because the transfers beyond them are neghgible. 
The Coulomb interaction term is considered in each vanadium site: 

'Hu ^ ^ Umm^ilma'^ilm—a ~l~ ^ ^ Umm'^ilma^ilm' —a ~l~ ^ ^ {^rnrn' Jrnrn') '^ilma'^ilm' c 
i,l,m,a i,l,mj^m,' ,cr 1,1,171^^171' ,<j 

— Jmm' (cli^^Ciim-aCiim'-a'^ilm'a — c\i^^c\i^_^Ciim' -trCUm' 1 (13) 

i,l,mj^m,' ,a 

where Umm and Umm' denote intra- and inter-orbital Coulomb repulsion, and Jmm' is the 
Hund's rule coupling. We have checked that the long-range interaction is irreverent for this 
system. The value of the interaction is shown in Table I. 

3. Hartree-Fock Calculation 

As a starting point of the PIRG calculation, we apply the Hartree-Fock approximation. 
The effective model is decoupled to a mean-field Hamiltonian in the usual manner. Concen- 
trating on homogeneous solutions, we take the order parameter as 

blmm'cr =^ y^Mllma'^i'-m'a) , (14) 

i 

where denotes number of unit cells. Assuming a specific SO as an initial state, we calculate 
the HF self-consistent equations until the energy converges. During the HF calculations, each 
SO is retained when the interactions are large, while the spin structure relaxes into a param- 
agnetic solution for small interactions. In Fig. 4, we show the relative HF energy for each SO 
measured from the ferromagnetic (FM) solution. It is consistent with the experimental result 
and the previous HF calculations^ that the G-type SO is the lowest in energy. The energy of 
the G-type SO solution is close to that of the C-type SO, which naturally explains the fact 
that YVO3 is located near the phase boundary between the G-type and G-type SO.^^ The 
HF ground state with G-type SO solution also gives the G-type 00. With these SO and 00, 
we also calculate the HF charge gap A^^ (Fig. 5), which gives 1.19eV for the realistic value 
(A = 1.0). Despite the HF method in general often overestimates the charge gap, our result is 
very close to the experimental one (~ 1 eV). Indeed, the gap amplitude is overestimated by 
2eV in ref. 9, although it also employs the Hartree-Fock approximation as well. The difference 
comes from the difference in the effective low-energy model. It is crucially important to derive 
a reliable low-energy model from the first principles. Our effective model obtained by the 
downfolding procedure has an advantage in reproducing both the order patterns and the gap 
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amplitude even at the HF level. We discuss the charge gap in detail below. 
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Fig. 4. Energy difference between the FM and various AF solutions as functions of A in the HF 
solution. The closed circles, diamonds and upper triangles represent the results of G— , C— and 
A— type AF solutions, respectively. The lower triangles is the results of the FM. 



4. PIRG Calculation 

Here we briefly summarize the procedure of the PIRG calculation.^ The ground state wave 
function is extracted from a proper initial wave function |<1>) by following the principle that 
limp^oo [exp [—rTiyf \^) generates the ground state where r represents a small number such 
that we can use the Suzuki- Trotter decomposition. The interaction terms are decoupled by the 
Hubbard-Stratonovich (HS) transformation. The path integral is expressed by the summation 
over the HS variables. The variational wave function is composed of L nonorthogonal Slater 
determinant basis functions as = ^i=i(^i\'Pi) ■ For a fixed L, the coefficients Cj and the 
choice of basis are optimized. Note that the case of L = 1 corresponds to the HF solution. 
The spatial and dynamical fluctuations are taken into account by increasing a number of 
Slater-determinant basis functions L. For sufficient large L, the ground-state energy decreases 
linearly as a function of the energy variance. This linear behavior ensures valid extrapolation 
to the full Hilbert space. 

For YVO3, the wave function of the G-type SO solution obtained by HF calculations as 
an initial state of the PIRG calculation gives the lowest energy state. We have performed 
calculations up to L = 192 oniV = 2x2x2 unit cells (32 vanadium atoms). We discuss 
system-size dependence later. In Fig. 6, we can see an example that this linear extrapolation 
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Fig. 5. The A dependence of the charge gap Ac for the HF {L = 1) solutions (circles) and the PIRG 
results (squares). The solid line is the least-square fit to the PIRG results. 



works well in a controllable way. We also check that the G-type SO and C-type 00 are 
preserved during the extrapolation process, while their amplitudes are slightly modified. The 
spin and orbital configurations in the ground state are illustrated in Fig. 7 

We calculate the total energy E{M-^,M^_) for the system with M = + M| electrons, 
where denotes number of electrons with spin a. Then we obtain the chemical potential as 
follows: 

E+ =E{4N + l,m + 1) (15) 
Eo =E{4N, m) (16) 

E_ =E{AN - 1, 4iV - 1) (17) 

AE E+ - Eq . , 

^+ = AM = 

AE Eo-E_ 
=AM = ' 

where Eq is the energy of the 1/3-filled system which corresponds to YVO3. The charge gap 
Ac is estimated as a difference between these two chemical potentials, 

Ac = Ai+-Ai_. (20) 

We note that the gap obtained by the above procedure is the indirect one that should give the 
lower bound of the direct gap. Figure 5 shows A dependence of the charge gap calculated by 
HF and PIRG. As the result of fluctuation, Ac of PIRG is reduced by 30-40% in comparison 
with the HF results. We also see that Ac behaves linearly as a function of A, which is similar 
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Fig. 6. The ground-state energy per unit as function of energy variance calculated by PIRG for (a) 
A=0.8; (b) A=1.0 and (c) A=1.2. Eq, denotes the extrapolated value of energy to the energy 
variance. 
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Fig. 7. (Color online): Ordered spin- and orbital- patterns in the ground state of the PIRG solution. 
The arrows represent the magnetic local moment at each vanadium atom. The orbital states are 
shown in the form of the spatial electron distribution. 

to the HF results. For the realistic value (A = 1.0), the charge gap is 0.70 it 0.07 eV, which 
seems to be slightly smaller than the experimental optical gap.^^ This is partly because that 
the PIRG calculation gives the indirect gap in general, while the optical gap is the direct one. 
The difference between the indirect and direct gap is estimated as ~ O.leV in the HF solution. 
Since the PIRG calculation is performed on a finite-size system, the charge gap may further 
be reduced with a finite-size scaling. Because of the limitation of computer resources, for the 
moment, results on finite-size scaling are not available. However, we have checked that the 
charge gap is insensitive to the system size in the HF calculations. It may suggest that the 
charge gap is determined from a local origin. We believe that the present result is close to the 
thermodynamic limit. Our result is quantitatively consistent with the experimental results as 
compared to those of HF(3.4eV)^ or GGA (0.009eV).^^ We stress that the present results are 
obtained from the first-principles calculation without any adjustable parameters. 

5. Summary 

In summary, the transition metal oxide YVO3 is investigated by the DFT-I-PIRG scheme. 
First we use the DFT method to construct the effective low-energy Hamiltonian which con- 
sists of the tight-binding part and the effective screened Coulomb interactions. The PIRG 
method is employed as a low-energy solver that can fully take into account spatial and dy- 
namical fiuctuations. The obtained spin and orbital state is consistent with the experiment. 
The indirect-charge gap is estimated to be 0.70 it 0.07 eV, which is smaller than the inferred 
experimental optical (direct) gap, but is consistent each other in terms of the difference be- 
tween direct and indirect gaps. It has prominently improved the estimation compared to the 
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LDA or GGA method. Our results are all consistent with the available experimental results. It 
is desired to measure the indirect gap to compare with the present result. The present appli- 
cation to YVO3 further proves that DFT+PIRG offers a quantitatively precise first-principles 
scheme for strongly correlated electron systems. The method establishes it as a standard 
scheme for reliable estimates of the electronic structure when the lattice structure is given. 
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Appendix: Transfer Integrals 

We present off-diagonal matrix elements in the transfer integrals of the tight-binding 
model. The indices i and j specify a unit cell and / and /' represent a vanadium atom in the 
unit cell. Between two vanadium atoms, the transfer matrix tjj;i'mm' has (m, m') elements, 
where m and m' denotes an orbital. Here we only show the transfer matrices with at least one 
element larger then O.OleV. Other transfer matrices are also included in the actual calculations. 

• i = j (intra-unit cell) 



t 



iillmm' 



/-0.27113 -0.00000 0.00000 \ 

0.00000 -0.21752 0.00000 

0.00000 -0.00000 -0.10540y 

/0.1085 



t 



iil2mm' 



t 



0.0407 
\0.0394 

/ 0.0366 
-0.0983 
\-0.1255 

'^0.0073 
= 0.0075 
y0.0288 



0.0407 


0.0392 \ 


-0.0362 


-0.0108 


-0.0110 


0.1323 y 


0.0144 


0.0031 \ 


-0.1206 


-0.0259 


0.0303 


-0.0217y 


0.0010 


-0.0044\ 


0.0027 


0.0067 


0.0047 


0.0082 J 



(A-1) 



(A-2) 



(A-3) 



(A-4) 



Tj = (a, 0,0) 
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t 



ijllmm' 



/-0.0156 -0.0058 -0.0077\ 

-0.0058 -0.0102 0.0146 
^-0.0078 0.0146 0.0042 j 



ijlSmm' ^ij24mm' ^iil'imm' 

/-0.0145 0.0148 0.0351 \ 

0.0108 0.0017 0.0003 
\ 0.0013 0.0022 -0.0058/ 



(0,6,0) 



ri - Vj = (a, b, 0) 



n - Tj = (a, 0, c) 



"^i - fj = 0' -c) 



t 



i+xj 14mm' 



^ijllmm' 



/-0.0052 0.0069 0.0002 \ 

0.0069 -0.0420 -0.0009 
^ 0.0002 -0.0009 0.0202 j 

^ijlSmm' ^ij24min' ^iillimm' 



/-0.0134 -0.0003 0.0149 \ 



tijll 



Ijllmm' 



-0.0003 0.0013 0.0040 
0.0149 0.0040 -0.0092 / 



^ijl2min' ^ij^imm' ^iil2mm' 
^ijlAmm' ^ij32mm' ^i-\-xjl4mm' 



^ijlSmm' ^ij2imm' ^iilUmm' 



ti 



ijlimm' ^jj23mm' ^i+xjUmm' 



t 



ijlimm' ^iiUmm' 



(A-5) 
(A-6) 
(A-7) 



(A-8) 

(A-9) 
(A-IO) 



(A-ll) 

(A-12) 
(A-13) 

(A-U) 
(A-15) 

(A-16) 



t 



ij23mm' — tiilimm' 



(A-17) 
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Full Paper 



(A.18) 



ri - Tj = (0, b, -c) 



''ij23mm' ''i+xjlimm' 



(A.19) 



(A-20) 



rj = {-a,-b,c) 



(A-21) 
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